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Abstract 

The present review is devoted to the problems of finite-size scaling due to the presence of long- 
range interaction decaying at large distance as l/r d+cr , a > 0. The attention is focused mainly 
on the renormalization group results in the framework of 0(n) c/? 4 - theory for systems with fully 
finite (block) geometry under periodic boundary conditions. Some bulk critical properties and 
Monte Carlo results also are reviewed. The role of the cutoff effects as well their relation with 
those originating from the long-range interaction is also discussed. Special attention is paid to 
the description of the adequate mathematical technique that allows to treat the long-range and 
short-range interactions on equal ground. The review closes with short discussion of some open 
problems. 
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I. INTRODUCTION 



Any system which has a finite size L in at least one space dimension we will call a 
finite-size system or a system with restricted geometry. In such systems singularities in 
the thermodynamic functions at the critical point may occur only in the thermodynamic 
(bulk) limit, taken in at least d\ dimensions (di is the lower critical dimension). Mainly three 
specific geometries (with periodic boundary condition imposed along the finite dimensions) 
are of particular interest: i) the fully finite cube L d , ii) the <i-dimensional layer (or film) 
L 1 x oo d_1 and iii) the infinitely long cylinder L d ~ l x oo. Further in this review we will not 
discuss the latter two cases. So far, in the context of long-range (LR) interaction they are 
studied only in the spherical n = oo limit. For the corresponding results one may consult 
the review [l| and the recent monograph |2|. 

The geometry of real systems and Monte Carlo calculations usually correspond to the 
fully finite size case. During the last two decades the study of systems with restricted 
geometry has undergone an extensive development and gains still growing importance for 
the theory of critical phenomena 0, H 0, 0, E|. Generally speaking the critical behavior 
depends essentially on geometry, boundary conditions, and on the universality class of the 
bulk system. 

It is well known that the universality class to which the critical behavior of a bulk system 
at a second order phase transition belongs depends upon the dimensionality of the space d, 
the number of the components of the order parameter n, the symmetry of the Hamiltonian 
(either in spin-space or in coordinate space) and the interaction potentials 0, H, 0, Eof . 
One of the most commonly studied interaction potentials, is the one corresponding to LR 
ferromagnetic interaction decaying algebraically with the spin separation r as r~ d ~° ' . The 
parameter a > controls the range of the interaction. The investigation of such systems 
was initiated by Joyce in his paper on the phase transition in the ferromagnetic spherical 
model^]. The interest in this type of interaction is tightly related to the exploration of 
the critical behavior of systems with restricted dimensionality in which no phase transition 
occurs otherwise (see 0, 0, 3| and references therein). Here the condition a > is 



needed to avoid an ill-defined thermodynamic limit. In the limit a — > + the interaction 
(after appropriate renormalization) equals the equivalent-neighbors interaction [15| . 

The opposite case, corresponding to a LR interaction with —d < a < 0, sometimes called 
also nonintegrable interaction, has been studied in connection with the so called "non- 
extensive thermodynamics" (see and references therein). In spite of the recent growing 
interest, the case of nonintegrable interaction is beyond the scope of the present study. 

The LR interaction, r~ d ~ a , enters the expressions of the theory only through its Fourier 
transform. The corresponding small q expansion of the Fourier transform has the general 
form 

v(q) = v + v 2 q 2 + v a q a + w(q) < a ± 2, (1) 

with w(q)/q a — > 0, for q — > 0, i.e in the long-wavelength approximation. Note that the case 
a = 2 leads to logarithmic factors in (JTJ which are not allowed to enter the analysis below. 
However, further on we will formally relate a = 2, to the short-range (SR) interaction since 
then is the Fourier transform of an interaction decaying exponentially with distance. 
Depending on whether a ^ 2 in (0) we will speak about leading or subleading LR term 
respectively. Since the SR term ^ q 2 in (P) always exists the problem in the momentum 
space must be considered with necessity as an interplay between SR and LR effects. 

This review is arranged as follows. In Section |H] we summarize the results on the critical 
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behavior of bulk systems with the LR interaction. Some general remarks on the finite-size 
scaling (FSS) for such systems are presented in Section IIHI We review some results on the 
finite-size critical behavior of systems with leading LR interaction in Section IIVI and for 
the subleading one in Section IV II Section |V| specially deals with Binder's cumulant as an 
utilized tool to compare the theory to data from Monte Carlo simulations. Some remarks on 
the cutoff effects are devised in Section IVllI In Section IV 1111 we deliberately emphasize some 
unsolved issues. In El we present some mathematical details in calculating lattice sums in 
the case of LR interaction. 



II. BASIC RESULTS IN THE BULK CASE 



The results concerning the bulk critical behavior in the simplest case of the spherical limit 
was analyzed in details in [TtJ (for a recent review, see Chapter 3 of These results were 
generalized to the 0(n) vector <£> 4 model by means of perturbation theory in combination 
with the renormalization group (RG) technique near the upper critical dimension 0, 0, 
El Ell d = 2a, lower critical dimension 0, El, d = a, and the l/n -expansion 



|27L l28j] . Computer simulations also contributed to the exploration of the critical properties 
of such systems @, H H H H H . 

In the case of LR interaction, some RG predictions, e.g. using the e-expansion, can be 
verified in an ideal testing ground since the value ofe = 2cr — d or e = d — a would be 
small enough for integer values of the dimensionality. In this context the outcome of the 
computer simulations, obtained by means of the Monte Carlo method, concerned mainly 
Ising systems (n = 1) with classical (mean-field) critical behavior (29I. l34|. i.e. < a < d/2 
with d = 1,2,3, and nonclassical critical behavior jlH, HH 35||, i.e. d/2 < a with d = 1, 2 . 



Some comparisons with rigorous results were obtained for low dimensional systems in |12j,[R 
and for d- dimensional systems in |36| . 

It is worth mentioning that the critical behavior depends strongly upon the interaction 
parameter a. When a < 2 the expansion has been used for detailed investigations of 
the critical behavior of 0(n) models including questions like the a, d and n dependence of 
the critical exponents and critical amplitude ratios, as well as for determining the universal 
scaling functions. In this case, the LR term is leading and the critical exponents of the 
system are a dependent. By increasing a, a crossover from LR critical behavior to SR one 
takes place. This issue has been the matter of a long-standing debate in the literature. As 
a measure for the crossover one may consider the Fisher exponent 77. 

First, it has been argued, that in the interval d/2 < a < 2 and above the critical 
temperature T c , the critical exponent 77 is equal to 2 — a with no corrections to order e 2 and 
e 3 (at least), see[3], and to 0(n~ 2 ), see |27], 0]. It is reasonable to believe that 77 "sticks" 
to this value to all order in e jjj. So 77 = 77(0") as a function of a, is not continuous at the 
crossover point" a = 2, since the SR exponent 77(2) 7^ 0. 

Later it was pointed out, in the cases of (4 — e) dimensions 0] and in (2 + e) dimensions 
(25l | that the crossover from LR to SR critical behavior is shifted and occurs at a " critical" 
value of a given by a c = 2 — t?(2). As a result rj = 77(c) is a continuous function in a at 
a = 2, since f](<J c ) = v(^)- Further support for this assertion was obtained in the framework 
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of different perturbation schemes in |20|, |21( and [23] . In terms of the RG language one can 
say that, by crossing the border point a c the corresponding LR and SR fixed points exchange 
stabilities and in a non vanishing range 2 — t?(2) < a < 2 LR perturbations are irrelevant (see 
also 26] where the problem in the more general case of (n, d, a) space was considered) . This 
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statement, which can be determined also from more general considerations (see for example 
p. 71 in [10i |) seems to be accepted in the literature. However, some doubts arise since it 
conflicts with the opposite statement of 36|| at least for the non Ising case n > 2 obtained 
on a most rigorous level. Furthermore there is the criticism in HI a gainst the use of (2 — a) 
as small perturbation parameter. If one reconsiders the problem with (2 a — d) as a small 
parameter in conjunction with the use of the perturbation scheme of 2^|, previous results 
that LR fixed point is stable up to a = 2, are restored. A recent attempt to reconsider 
this issue using Monte Carlo simulations |35( shows unambiguously, a crossover at a c , in 
agreement with 0, H^, but there only the Ising case (n = 1) is considered. Consequently, 
the situation seems to be settled only for the case n = 1. 

When a > 2 one usually considers the model as equivalent to a model with SR interaction, 
since it is widely accepted [HI |25l |2(| |2?J that, the LR terms (i.e. the third and etc. in 
Eq.((TJ)), do not contribute to the critical behavior of the bulk system, consequently this 
term has been always omitted in the computational analysis. Indeed, in this case, the 
critical exponents does not depend on the parameter a. 



III. SOME GENERAL REMARKS ON THE FINITE-SIZE SCALING 

Scaling is a central idea in critical phenomena near a continuous phase transition and 
in the field theory. In both cases the singular behavior emerges from the overwhelming 
large number of degrees of freedom, corresponding to the original cutoff scale, which need 
to be integrated out leaving behind long wavelength with smoothly varying momenta. This 
behavior is controlled by a dynamically generated length scale: the bulk correlation length 
£. Such a fundamental idea is difficult to test theoretically because it requires the study 
of a huge number of interacting degrees of freedom. Experimentally, however, one hopes 
to be able to study scaling in finite systems near a second order phase transition. Namely 
the system is confined to a finite geometry and the FSS theory is expected to describe the 
behavior of the system near the bulk critical temperature (for a review on the FSS theory 
see 0, H In a few words the standard FSS is usually formulated in terms of only one 
reference length - the bulk correlation length £. For a system with finite linear size L, the 
main statements of the theory are: 

i) The only relevant variable in terms of which the properties of the finite system depend 
in the neighborhood of the bulk critical temperature T c is L/C,. 

ii) The rounding of the the thermodynamic function exhibiting singularities at the bulk 
phase transition in a given finite system sets in when L/C = 0(1). 

The tacit assumption is that all other reference lengths are irrelevant and will lead only 
to corrections towards the above picture. Moreover, the crucial point in the FSS theory is 
that we always assume: the finite size L of systems under consideration and the correlation 
length C are large in the microscopic scale. This means L ^> a and C 3> a, where a is 
the lattice spacing. The scaling limit of a quantity is called its value when all corrections 
involving the ratios a/L and a/C are neglected. One postulates FSS limit as: 

C L L 

> oo, > oo, — = const. (2) 

a a 4 

In the literature there exist two conceptually different approaches concerning the scaling 
regime. One inspired from the condensed matter theory uses lattice language and models, in 
which a is fixed (while L and C are large, going to oo). The other one is inspired from particle 
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physics. It uses continuum field theory models, when L and £ are fixed (while a is infinitely 
small, going to zero). For studying a dimensionless quantity when all the dependence on 
a, L and £ is through the ratios a/L and a/£ the two approaches seem to be equivalent. 

We shall be interested in the case where only the bulk system displays a phase transition. 
In such system there is no fixed point for finite L and so no crossover to any other fixed point 
as a function of L/£ can takes place. This means that an L- independent RG procedure can 
be used. In the case of LR interaction there are additional reasons (related with possibilities 
of phase transitions in lowered dimensions) to study the case where finite system as well 
the bulk displays a phase transition. However this more sophisticated case has been studied 
only in the SR case |37l |. 

So, we review the FSS properties of fully finite systems in the framework of the continuous 
field theory. The effects of the cutoff A = n/a related with a fixed will be discussed in some 
details below. 

In case of LR interaction there are two aspects of the theory which should be mentioned 
here. First, in addition to the SR term |q| 2 in the propagators of the Feynman diagrams 
we have also the LR term |q|°". The last term causes peculiar mathematical problems 
jlH 3^, 4(| concerning the evaluation of the lattice sums over q. Among the various methods 
of accomplishing these summation some details of the method we prefer is presented in 1X1 
where its advantages are also discussed. The second aspect arises specifically because the 
basic quantity of the theory £ is unambiguously defined for SR interactions. For the LR case 
it has to be defined in a more appropriate way jZJ]. The standard definitions of the bulk 
correlation length are based on the asymptotic behavior of the pair correlation function at 
large distance. The two most commonly used definitions of the bulk correlation length are 
unambiguous in the case of exponential decay of the bulk pair correlation function Goo(R; T) 
with the distance R = |R|. One defines 

£i(T) = - lim[i2/]nGoo(R;T)]. (3) 

Alternatively, one may consider the second moment of the bulk pair correlation function 
and define the effective correlation radius, 
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Y / R 2 Goo('R;T)/G 00 (R ] T) 



R 



1/2 



(4) 



If the pair correlation function decays exponentially with distance, the two definitions are 
equivalent. A distinctive feature of the LR interaction is that the function G^R; T) decays 
as R~ d+cr when R — > oo, for T > T c (see e.g. jrH]) and both definitions yield £i(T) = £2 CO = 
oo. One can overcome this difficulty by following the approach proposed in [l| (see also p. 138 
of 0) where instead of £ some bulk characteristic length A(T) is used. This characteristic 
length determines the length-scale of variation of the correlation function and diverges at 
the critical point. 

As in the bulk limit 0(n) models are the most often used as laboratory tools on the basis 
of which one studies the scaling properties of finite-size systems. The exactly investigated 
cases are those corresponding to n = 1, i.e. for the Ising model (particular attention has 
been turned to the two dimensional case), and the limit n = oo, which includes the spherical 
model (see e.g. Chapter 5 in 0] and p. 141 in 0]). It is commonly assumed that the last 
one is the only model which combines exact solubility even in the presence of a magnetic 
field and direct relevance to the physical reality. For that reason it is especially suitable 
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for testing the FSS hypotheses. For n ^ 1, oo there are no exact results and the commonly 
preferable analytical method for the derivation of the properties of the corresponding models 
(like XY, i.e. n = 2, and Heisenberg, i.e. n = 3) is that of the RG theory. In addition 
an important amount of information for such systems is derived by numerical simulations, 
normally via Monte Carlo methods. 

A systematic and controlled field-theoretical approach to the quantitative computation 
of the thermodynamic moments was proposed in the middle of 8Q'sj42L 0] for studying 
FSS in the case of SR interaction. It is based on the idea of using a mode expansion, i.e. 
one treats the zero mode of the order parameter, which is equivalent to the magnetization, 
separately from the higher modes. The nonzero modes are traced over to yield an effective 
Hamiltonian for the zero mode. This method is used in combination with the loop expansion 
in the framework of the minimal subtraction scheme. The method is quite general and was 
proven to apply to a large extent to the investigation of FSS in systems with LR interaction 
and with finite number of components n of the spin vector 1441. However such systems have 
been the object of analytical investigations mainly recently |32l 45l . 46l 47| . For example, it 
is possible to perform the quantitative computation of the thermodynamic moments, usually 
used in numerical analysis. These moments are related to the Binder's cumulant B and to 
various thermodynamic functions like the susceptibility. 

The present wisdom is that in the finite size case the results are quite different depending 
on whether the LR interaction is leading or subleading. So, below we will consider separately 
the cases of leading and subleading LR interaction. 



IV. LEADING LR INTERACTION 

Recall that the critical exponents depend upon the parameter o < 2 controlling the 
interaction range. The FSS properties (under periodic boundary conditions) in the spherical 
limit are well established. It has been found that the scaling properties of the system with 
SR interaction remains valid also in the case of LR interaction (for a review see e.g Chapters 
4 and 5 in J2| and references therein) . For finite n a limited number of recent numerical 
results 3(1 22, Hi| , as well as few analytical works 3^, 3 3, Hil |46| became available. 



The nonlocal character of the LR interaction has been an obstacle in investigating the 
critical behavior by means of numerical methods. For that purpose different algorithms have 
been developed. Using the Ewald method to evaluate the energy of a given configuration, the 
critical behavior of such systems for different values of d, n and a has been investigated 3^ ■ 



Another approach, analogous to FSS, where the range of the interaction is cutoff at a certain 
value, whence the name "finite-range scaling", is also used.jjyj This method has been tested 
essentially on the one dimensional Ising model in the classical (mean-field) as well as in the 
nonclasical critical regimes. Results of the aforementioned approaches have been judged 
comparable to the theoretical predictions. These treatments required special efforts that 
restrict the consideration to very small systems. Recently, the problem was resolved by the 
use of cumulative probabilities within the Wulff cluster algorithm [29|, |34|. This method 
has been applied mainly to Ising systems in the mean-field regime. The FSS hypotheses 
have been tested in the nonclassical regime [3^, 34||. The common conclusion of all these 



approaches is that in the case d = a such systems exhibit a Kosterlitz-Thouless-like phase 
transition. 

Analytically the nonclassical case has been studied first for n = 1 3^ and for n ^ 1 



45l |46[. It has been found that, as for the bulk systems, the critical behavior depends on 
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the small parameter e = 2a — d. The scaling properties of the finite system are not altered 
by the presence of the LR interaction and can be written as 

X = L-^ /v T{tL 1/v ), 

where j x is the critical exponent of the observable X and t = (T — T c )/T c . Some thermo- 
dynamic quantities, like for instance, the shift of the critical temperature due to the finite 
size effects, the susceptibility and the Binder's cumulant at the bulk critical temperature T c 
'32I 45L and above it, ji^, 4f| as a function of y/e have been obtained. The d dependence 



of the finite size properties has been considered in details in reference [46] by means of a 
method, which consists of the use of the minimal subtraction scheme applied to a fixed space 
dimensionality (48[. 

It has been shown that the critical behavior of the system is dominated by its bulk critical 
behavior away from the critical domain and that the FSS is relevant in the vicinity of the 
critical point. A distinctive feature of the LR case is that for tL x l v ^> 1 the finite-size 
corrections are not exponentially small as in the SR case; they vary instead as power-in-L 
law, in both the spherical n = 00 0> 41] and n ^ 00 cases |45l E^. 



In the behavior of the coupling constants during the crossover from LR to SR in- 
teractions i.e. the limit a — ► 2, was considered. To one-loop order it has been concluded 
that the renormalized values of the temperature and the coupling constant are continuous 
functions of a. 



V. BINDER'S CUMULANT 



As it was mentioned above a quantity of central interest in the study of the critical be- 
havior of finite size systems is the Binder's cumulant ratio B. Here we present the analytical 
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result for B (B = 1 — M 4 /3Mf , where M 2n is the n-th moment; for some details, see 42] and 



also [32j,|45j). Close to the critical point we have 



B 



1 - 



n 



r 2 Q») 



12 (J(n + 2)) 



r(> + 6)) r(l(n + 2)) _ r(|(n + 4)) 
r(I(n + 4)) + r(|n) 



rQ(n + 6))r(l(n + 2)) r 2 (i(n + 4)) 
r(J(n + 4))r(|n) P(J(n + 2)j 

At the fixed point z is given by 45] (e = 2a — d) 



n 




n + 8 



n 



n + 8 



hxy +2 a - 1 e 



„_! n + 2 



n + 8 



T{a)F 2a ^ (y) 



e2°->yT{a)F^ a 



Here, we have introduced the scaling variable y = tL l l u and the function 

„<t/2 



F d ,a (y) 



dxx 2 1 Es. 



yx 



(2tt) ct 



6 d (x/n) 



ir\ d / 2 
xJ 



(5) 



(6) 



(7) 



In Eq.© E ai p(x) and 9(x/ir) are the generalized Mittag-Leffler and the reduced Jacobi 63 
functions, respectively (see also 1X1 for details). The function F d>c (y) is well known in the 
theory of FSS. 
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FIG. 1: The cumulant ratio B as a function of the interaction range a for d = 1. The case 
< cr < d/2 corresponds to classical (mean-field) regime. The Monte Carlo results "•" follow from 
for < a < 0.5 and from 32] for 0.5 < a < 1. 

Finally, let us notice that one can easily see that the expression for z as a function 
of y verifies the FSS hypotheses and, consequently, all the thermodynamic functions, which 
indeed are z dependent, do (Hf. At the critical temperature T c (i.e. t = 0, and so y — 0), 
we obtain 

r n + 2 



r< 



F 2ff)ff (0)+O(e) 



(8) 



where the coefficient -F 2o - jCr (0) appearing in the right hand side of (jHJ) can be evaluated |4£ 
analytically as well as numerically for different values of the interparticle interaction range 

r 2c (1/2) , 



F 2 a,a (0) 



4C (1/2)/? (1/2) 
-4.82271993, 
-8 In 2, 



a = 1/2, 

<7=1, 

a = 3/2, 
a = 2. 



(9) 



Here C(x) is the Riemann zeta function with ( (|) = —1.460354508... and (3{x) is the analytic 
continuation of the Dirichlet series: 



oo 



-r 



with /5 (|) = 0.667691457.... Note that the function F 2u , a (0) increases as the parameter a 
vanishes. 

Numerical values for the Binder's cumulant ratio Q can be obtained by replacing the 
value of Zq form (jHJ) and taking some specific values of the small parameter e. The behavior 
of the universal constant B for the cases d — 1, 2 is presented in Figures [T] and |2l respectively. 
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FIG. 2: The same as Fig. HJfor d = 2 



Note that the scaling variable z is proportional to y/e as it was found previously (see 4^ 
for example) in the case of SR forces. Furthermore for n — 1 it is in full agrement with the 
result obtained in reference devoted to the exploration of FSS in 0(n) systems with LR 
interaction. Notice that in |32j| (see also jsij) the mathematically ensuing pertinent integrals 
have to be evaluated only numerically, due to the choice of a parametrization, e.g. (|A9J) . 
that does not reduce the <i-dimensional problem to the effective one- dimensional one. The 
approach proposed in [45| based on the parametrization (jAlOj) is more efficient in the sense 
that the corresponding final expressions can be handled by analytical means |49| . 

It has been found 32], using the Monte Carlo method, that the amplitude ratio Q = 
M\jM{ (which is related to B) is a linear function of e (see also Fig.l and Fig.2), while the 



analytical result (jHJ) shows an expansion in powers of y/e |22|, |4j|. It is possible that higher 
orders in e could improve the result. Better agreement could also be obtained by using the 
method developed in Ef 



VI. SUBLEADING LR INTERACTION 



Let us first note that the most prominent example for such kind of interaction in the case 
d — 3 is the van der Waals interaction 1/r 6 for which we have a = 3. One can easily see 
that in this case the Fourier transform is indeed of the type (jT|. 

The fact that the critical behavior is not affected by the LR interaction if a > 2 is well 
established in the bulk case. Analytically, in the framework of n = oo model it has been 
shown j^j that such a statement is incorrect for finite-size systems because of finite-size 
contributions due to the subleading, a > 2, term in the interaction. It has been shown also 
that the same remains true for a finite number of component by means of the RG techniques 
j4?|. The authors of j4?J have investigated the FSS behavior of a fully finite 0{n) system 
with periodic boundary conditions and in the presence of a LR interaction that does not alter 
the SR exponents of its critical behavior. The small |q| expansion of the Fourier transform of 
the interaction v (q) is supposed to be of the form (JTJ) with 2 < a < 4. For such a system, it 
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has been demonstrated that all the thermodynamic functions can be expressed in a scaling 
form as 

X = L-^ v F{tL l l v ,bL 2 - a -v), 

where b is a model (nonuniversal) constant. Note that one needs two scaling variables in 
order to describe in a proper way the finite size behavior of these quantities. 
At the critical point the Binder's cumulant is given [4?} by (jSJ) with 

z = -V^—^L [In 2 + bzu a L^ + 0(e)] , (10) 

7T a/72 + 8 

where w a = (2n) a - 2 (l - 4 <T / 2_1 )C(1 - cr/2)C(2 - a/2). So, in this case B is not an universal 
constant. The behavior of B as a function of system size L is presented in Fig. El Equation 
(TTTJ|) is a generalization for the case under consideration of the result obtained in ji^l for SR 
interaction. 

When tL 1 ^ > 1, it was found that the susceptibility approaches its bulk value not in an 
exponential-in-L, as it was commonly believed to be the case for systems with short-range 
critical exponents, but in a power-in- L way of the order of bL~( d+a \ The last goes beyond 
the standard formulation of the finite-size scaling, but is completely consistent with the 
intrinsic large-distance power-law behavior of the correlations in systems with subleading 
LR interactions. In the spherical limit such large-distance power-law behavior is shown 
exactly in [5l| . 

For an arbitrary value of <r, away from the critical point, i.e. in the region tL l ' v > 1, 
the Binder's cumulant ratio is given by B = 1 - |(1 + |), with finite size correction falling; 
off in a power law. This result corresponds to a n-dimensional Gaussian distribution for n 
independent components of the vector variable. Obviously, all the values lie in the interval 
from B = (Ising model, n — 1) to B = 2/3 (spherical model, n = oo). 
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VII. LATTICE AND CONTINUUM MODELS AND THE EFFECT OF CUTOFF 



In the above, we have concentrated our attention on how to separate the size dependence 
in the continuum (scaling) limit, where the linear size as well as the correlation length tend 
to infinity, but keeping their ratio a finite quantity. In this theory the cutoff is sent to 
infinity and the lattice spacing completely disappears. This is precisely the regime in which 
one expects FSS to hold |42|,|43[ and this statement, as it was pointed out, does not depend 
on the range of the interaction. 

In order to clarify the effect of a (sharp) cutoff it is sufficient to consider the pure SR 



case. The finite cutoff effects violate the FSS, see |52|,|53[ and references therein. The result 
is a system with finite-size behavior with leading nonscaling term going as an inverse power 
law in L that depends also on LA. In [H^l this is shown to be tightly related to lattice 
effects in the system in conjunction with the long wavelength approximation. Indeed this is 
not in conflict with the arguments proposed in |42l . |43| where the limit of infinite cutoff is 
considered. Later, in one argues : i) the violation of scaling is related with the second 
term in (jA4j) which in the case a ^ can not be omitted in the region L/£ ^> 1, ii) the 
effect of the sharp cutoff is similar to that found if one considers the effect of subleading LR 
interaction. Obviously, this effective LR nature of the finite-size behavior should generate 
difficulties in the analysis of the Monte Carlo data from simulations not only for SR, but 
for LR systems as well. The last follows directly from (jA13|) . Recently, the reconsideration 
of the problem in has shown that the source of these obstacles is the artificial cusp-like 
singularity 



dv(q) 



dqi 



_^ dv(l) 



dqi 



d (11) 



at the border of the Brillouin zone (A ~ 1/a). As a result, if one treats properly the effects 
generated by the momentum cutoff in the Fourier transform of the interaction potential both, 
the lattice and continuum models may produce results in mutual agreement independent of 
the cutoff scheme Ifk 



VIII. UNSOLVED PROBLEMS 



First, an area where the bulk RG predictions HI HI must be tested by means of numerical 



methods is the case of n > 1. The strong indicatio n I36fl that exactly here some problems 
exist, remains still valid. The numerical results of [35j] shed some light on the Ising case 
n = 1 but do not solve the longstanding controversy about the boundary between LR and 
SR critical behaviors function of a. 

In this review the FSS properties of systems with no crossover from the bulk fixed point 
to any other as a function of L/£ are considered. Let us note that within the SR models this 
issue is studied in j3?J ■ The consideration of this crossover in the LR is missing. On the other 
hand, e.g. systems with slab and cylinder geometries and n ^ oo have been studied only in 
the SR case. The first one is related to the statistical-mechanical Casimir effect in fluctuating 
systems |55| and recently there has been an upsurge in theoretical investigations 0, IE H3 5 E3- 
The second one deals with the so called transfer matrix (or Hamiltonian) formulation of the 
problem (see also Chapter 36.4 of [9]). A study of both cases in the context of the LR 
interactions is a provocation for the theory and is of undoubted experimental interest. 

Some problems exist with the comparison of the analytical results and the numerical 
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data. It has been shown [321 (see also |3^|) using Monte Carlo simulation (for n = 1), 
that the Binder's cumulant ratio is linear in e. The analytical evaluation (jSJ) for the 0(n) 
symmetric <y9 4 model, however, showed that it is linear in ^ft. A possible way to resolve this 
controversy between the Monte Carlo method and the analytical results is to carry on finite- 
size calculations to higher loop order j^. This could ameliorate the analytical results, which 
would be comparable to those obtained by numerical simulations. However we would like to 
mention that higher loop corrections that are dealt with through the minimal subtraction 
scheme and the e-expansion are not done, even for the more simple case of SR interaction. We 
have witnessed that two-loop calculations can be applied in many other investigations like 
the problems taking into account disorder effects in finite size systems jH^. Here, if n = 1, 
the one-loop fixed point is degenerate and the first interesting results can be obtained only 
if one considers two loops. 

As it is discussed in the previous section in the long wavelength approximation if 
one goes beyond the continuum field theory, the lattice effects in conjunction with the sharp 
cutoff generate nonuniversal finite-size terms. The suggestion that such effects are artificial 
in the case of SR interaction [54| seems to be true also in the LR case. How to avoid the 
problem in the framework of a concrete cutoff scheme in the case of LR interaction is still 
an open problem. 

In recent years considerable attention has been paid to the critical dynamics of systems 
with LR interaction It would be interesting to extend the corresponding results to the 
finite-size case. Let us note that for the case of the SR interaction the theory of finite size 
effects in critical dynamics was developed in [5^(see also Chapter 36.6 of 0). 

The present review is devoted to the classical critical phenomena. However, another 
interesting field is closely related to the extensively investigated field of quantum critical 
points, i.e. phase transitions occurring at zero-temperature jH^]. Let us just recall that 
in systems showing quantum critical behavior the temperature plays two different roles. 
For temperatures low enough, quantum effects are essential. In this case the temperature 
affects the geometry to which the system is confined adding a "new" size to the Eucludean 
space-time coordinate system. By raising the temperature, the system is driven away from 
the quantum criticality. At high temperatures, however, the size in the "imaginary-time" 
direction becomes irrelevant in comparison with all length scales in the system. In this case 
we have a classical system in d dimensions and the temperature is just a coupling constant 
in the classical critical behavior. In this context we find it useful to investigate the quantum 



models, with LR interaction, e.g. considered in |60|, |61| in the large n limit or in |62| in the 
mean field and one-loop RG theory. 



APPENDIX A: THE ORIGIN OF THE MATHEMATICAL DIFFICULTIES IN 
EVALUATING THE LATTICE SUMS 

In order to introduce the reader to the problems let us first consider in some details the 
SR case on a d- dimensional hypercubic lattice Z d with N = Nq sites and periodic boundary 
conditions. The sites are given by 

r = a(nx, . . . ,n d ), (Al) 
where rii ranges over all distinct integer values modN , for i = 1, . . . , d. 



12 



Let us consider the dimensionless expression , entering most of the mathematical analysis 
of the scaling properties of finite systems 



G (N ,a\d) 



(amN ) d J 



m 



_ d No/2-1 

d 



or — — m + Q 

nx=-N /2 n d =-N /2 |M| 



(oiV, 

30 

exp(— t) 



No/2-1 

E 




/ 2im 
V amNo 



dt, 



(A2) 



where Oj = and N is even. For the sum 



QN (t) 



, l v Wo/2-1 

(amiVo) exp 

u/ n=-N /2 



( 2 



7rn 



amN 



(A3) 



using the result of [5Cj, we have 
QN {t) - 



+ 



/47rt 
1 



erf 



7Tt 1 /2 



am 



27T 2 t 1 



3 am 
^exp[-Z 2 (amiVo)74t] 



exp 



-TlH 



1 

am 



(A4) 



valid in the large N asymptotic regime. The first and the second terms in the above equation 
are size independent and are precisely the infinite volume limit of Qjv (t). The ultraviolet 
divergences which may appear in the theory when the lattice spacing vanishes are related 
with the first term, for which one must perform the required subtractions (see, e.g. p. 177 
of 0). The second term depends on the ratio £/a = I /ma and so if a — > or £ — > oo it is 
of order 



O (£/a)te 



-nHtf/af 



(A5) 



In the continuum limit such terms are exponentially small and must be omitted. 

The different finite-size regimes are governed by the ratio £/L in the third term. They 
are independent of the microscopic details, e.g. lattice spacing a. When £ 3> a, the uni- 
versal properties can be described by a continuous field theory. Let us recall that the finite 
linear dimension L = Noa, in the case of continuous finite volume means that a — > and 
simultaneously N — > oo. 

In the continuum limit, since erf (nt 1 ^ 2 / am) = 1, we end up with the result 



G(N , 0|d) = G(oo, 0\d)\ L=oo + g{L\d), 
where the size dependence is contained in the function 



(A6) 



g(L\d) 



1 



^d/2 
1 



exp(— t) 



t d/2 



exp(-m 2 \\(d)\ 2 L 2 /4t)dt 

\(d)=i 
exp(— t) 



(4vr) rf / 2 J t d / 2 







m 2 L 2 



Ant 



- 1 > dt, 



(A7) 
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and where 



-rtxn 



(A8) 



The bulk part G(oo, 0|g?)|l =oo contains poles at d = 2,4, The finite size correction g(L\d) 

combined with the corresponding contribution from the counterterm yields 4^, 43[ (to one 
- loop order) a finite-size shift of the "mass term" m 2 . 

In order to investigate the FSS properties of systems with LR interaction, one can use a 
suitable mathematical method allowing to simplify the analytical calculations. In the above 
case of SR interaction it has been possible to replace the summand in eq. (jA2|) by its Laplace 
transform 



E 



1 



m 2 + q 



dtexp(— m 2 t) 



J^exp(-g 2 t) 



(A9) 



where q and q are d-dimensional and one- dimensional discrete vectors, respectively. This is 
the so called Schwinger parametric representation. The aim of this replacement is two fold: 

i) to reduce the d-dimensional sum to the corresponding effective one- dimensional one, and 

ii) to give to the dimensionality d the status of a continuous variable. In the case with the 
a term leading in (JT|. one cannot just use the Schwinger representation in its familiar form. 
In[39| the following generalization of (|A9|) has been suggested 



E 



1 



m 2 + |q| CT 



m 



dtQ CT {m A/a t) 



E 



exp(— q 2 t) 



(AlO) 



where the function Q a {t) for < a < 2 is given by 



and 



Qa(t) 



QM 



dyexp(-ty)Q a (y), 



sm(cnr/2)y' 



a/2 



it 1 + 2y a / 2 cos(cr7r/2) + y° 



From (|A10|) and (|A11|) for the summand in (|A10|) one obtains 



dtQ a 



m 



tm 2 + |q| 2 



(All) 



(A12) 



(A13) 



The integral representation (|A13j) illustrates the mathematical difficulty appearing in the 
LR case. It is shown that by an additional integration the problem can be effectively reduced 
to the SR casej39|. The behavior of Q c {x) as a function of a (see Fig. @| shows another 
aspect of the nonlocal character of the LR interaction l/r d+cr . 

The function Q a {x) is related [63| to the entire function of the Mittag-Lefner type defined 
by the power series 

(A14) 



E, 



a ,f3 



fc=0 



T(ak + [3) 



a > 0. 



For a more recent review on these functions see 



It has been demonstrated 



Q lT (x)=x^ 2 - 1 E aM2 (-x^ 2 ) 



that 
(A15) 
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a = 2.00 
a= 1.90 
a = 1.75 
a = 1.50 



3 - 



2 - 



0.5 



1.5 



FIG. 4: The dependence of kernel (|A12D as a function of x for different interparticle interaction 
range a = 1.50; 1.75; 1.90 and 2, 00. 



and all the specific features of the finite size properties of systems with LR interaction are 
a result of the analytical properties of the Mittag-Leffler type functions (for some examples, 
see e.g. Chapters 5 and 6 of 0). 

In the case with both SR and LR terms the corresponding expressions are significantly 
more complicated. When the LR term is of the peculiar type a = 2 — la (and the parameter 
a — > + ) considered in the context of the bulk crossover the following replacement of the 
SR propagator takes place HF 



1 



m? + |q| 2 



E 



(-v*y\ 


q| 


2Z(l-a) 


(m 2 + 


|q 


|2)1+Z 



(A16) 



In order to reduce the problem of evaluating the asymptotic behavior of the sum over q 
to the corresponding one- dimensional sum in the right hand side of (|A16[) . the following 
identity has been used 3 





q 


21(1- 


-a) 


(m 


2 + 


q 


2U+Z 



T(l + Za) 



dxx 



la 



iFi(l + I; 1 + la; -m 2 ) exp(-|q| V). (A17) 



Here, \F\ is the degenerate hypergeometric function (Kumer's function). The same iden- 
tity (jA17|) has been used in |50j to analyze the finite-size behavior of a propagator with SR 
and subleading, i.e. o > 2, LR interaction (treated as a perturbation) where terms similar to 
the r.h.s. of (jA16|) also appear. In the last case an expression for the entire propagator (i.e. 
SR and LR interactions treated on equal ground) can be obtained via contour integration 
on the complex plane |Ho|. 

The identities (jAlOj) and (jA17|) demonstrate that LR case can be effectively reduced to 
the SR case with an integration over an additional parameter. So, all the mathematical 
machinery developed for the SR case may be used without further complications. 
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